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When a cold atomic gas is illuminated by a quasi- resonant laser beam, light-induced 
dipole-dipole correlations make the scattering of light a cooperative process. Once a fluid 
description is adopted for the atoms, many scattering properties are captured by the 
definition of a complex refractive index. The solution of the scattering problem is here 
presented for spherical atomic clouds of arbitrary density profiles, such as parabolic densi- 
ties characteristic of ultra-cold clouds. A new solution for clouds with infinite boundaries 
is derived, that is particularly useful for the Gaussian densities of thermal atomic clouds. 
The presence of Mie resonances, a signature of the cloud acting as a cavity for the light, 
is discussed. These resonances leave their fingerprint in various observables such as the 
scattered intensity or in the radiation pressure force, and can be observed by tuning the 
frequency of the incident laser field or the atom number. 

PACS codes: 



1. Introduction 

Ultracold and quantum plasmas have shown an increasing interest for t heir analogies 
with cold atomic systems ( Mendonca fc Terc as 2013; Mendonca et al. 2008). In the dilute 
regime, direct particle interactions can often be neglected and the particles are interacting 
only via their common radiation field. Furthermore, cold atomic systems are often at the 
borderline between classical and quantum realms, allowing to investigate in a controlled 
way how cooperativity changes when a transition to a quantum regime occurs. A simple 
example is provided by a col d atomic cloud released from a magneto- optical trap and 
illuminated by a laser beam ( Bienaime et al. 20101 : Bender et al. 2010f) . Such a system 
manifests many interesting effec ts when the atomic cloud loses its granularity and can be 
described as a continuous fluid ( Prasad &; Glauber 201fj| ). As size, shape and density of 
the scatterer vary, a variety of radiating pattern s emerge, well known since the pioneering 
studies of Mie on extended dielectric particles (|Midll908t Ivan de Hulstlfl98lh . 

Although Mie scattering from an uniform dielectric sphere already represents a rather 
complex problem, light scattering by c old atoms adds new features that deserve dedicated 
studies ( Bachelard et aLlfeolll I20126I ). Firstly, changing the laser wavelength allows to 
tune the light-matter interaction: this allowed to observe experimentally the transition 
from single-atom s cattering to Mie scattering by a macrosco pic object, when the collec- 
tive effects appear ( Courteille et al. 201dL Bender et al. 2010h . Secondly, the weakness of 
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direct inter-particle interactions allow to probe the forces that light exerts on each indi- 
vidual atom: in particular, for ultra-cold atomic clouds the pattern of atomic recoil was 
shown to contain th e history of the interact ion, and to exhibit signatures of anisotropic of 
the Mie scattering (jBachelard et aZ.ll2012al ). Finally, various geometries, atomic densities 
and ordered or disordered structures can be generated by designing appropriate magnetic 
and optical traps for the atoms. This illustrates the flexibility of cold atom experiments 
to model and investigate phenomena from various fields, such as condensed matter and 
possibly plasmas. 

In this paper we present the solution of the Mie scalar scattering problem for spherical 
inhomogeneous atomic clouds, with finite and infinite boundary conditions. Beyond the 
classical Mie solution for hard and homogeneous dielectric spheres, it provides a realistic 
description of light scattering by dilute atomic clouds. We discuss the possibility for the 
cloud to act as a resonant cavity for the light, and the signatures of this behavior. 



2. The scattering equation 

Let us consider a sample of N cold two- level atoms with resonant frequency CJ a , 
linewidth V = d 2 Lj^/2Trheoc 3 and electric dipole transition matrix element d. When illumi- 
nated by an uniform laser beam with electric field amplitude Eq, frequency wo, and wave 
vector ko = (uJo/c)e Z) t heir response is described for weak incident fields by the follow- 



ing c oupled equations ( Scully et aZ.ll2006c ISvidzinskv et aZ.1 120081 120101 ICourteille et al 
2010h : 



N 



exp(ifc |r J - r m |) 
ikolr, - r m \ 



Pm{t) 



(2.1) 



with Ao = loq — co a , Tj the position of atom j, and /3j its complex probability amplitude 
to be excited at time t. The scattering kernel 



G(r,r') 



exp(ifco|r — r'| 
ik Q \r - r'l 



(2.2) 



has a real component that describes the collective (superradiant) atomic decay (jSvidzinskv et al. 
2008) , and an imaginary component that contains the collective Lamb shift due to light 



induced short-range interaction between the atoms (Friedbcrg et al ll973tlScullv fc Svidzinskvl 
20091 l2010h . The latter becomes significant when the number of atoms in a cubic optical 
wavelength, pA 3 , is larger than unity. Eq. (pOT) has been obtained assuming at most one 



atom is excited (jScullv et aD l200dlCourteille et aZjl2010h . However, this model has also 



been shown to describe the dynamics of coupled classical linear oscillators (jSvidzinskv et al 



|2010() , so that the many-body features of cooperative scattering by two- level atoms may 
be understood as a classical effect. Remark that this analogy holds only in the weak ex- 
citation limit of the atomic ensemble. Furthermor e, in this approach sho rt-range dipole 
interactions and polarization effects are neglected (|Friedberg et al. I ll973h . 

Neglecting granularity and isolating the self-decaying term —(T/2)j3j, E q. (12.11) takes 



the for m of an integro-diffcrential fluid equation for the complex field /3(r, t) (jBachelard et al 
201261) : 



dp(r,t) 
dt 



= iA 



2ft 



drV(r')G(r,r')/3(r',t) (2.3) 
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where p(r) is the atomic density. At steady-state, Eq. f|2 . 3[) yields 

drV(rO eXP ^f° |r ~| r/|) ^(rO = -(26 + i) /3(r) + (2.4) 
k \r - r'\ 

where we have introduced the normalized detuning S = Aq/T and excitation field /3(r) = 
(hr/dEo)f3(r). Let us remark that the kernel of the integral of Eq. (12.41) is the Green 
function for the Helmholtz equation: 

V*****?^--**-*) (2.5) 

and that (V 2 + k%) exp(ik • r) = 0. Then, applying (V 2 + fc 2 ) to Eq. (j2~4l) we obtain 
that /3(r) satisfies the Helmholtz equation 

[V 2 + *gmg(r)] /3(r) = (2.6) 

where mo(r) is the local refraction index of the atomic cloud: 

<«) 

Hence, the field /3(r) propagates as a wave in the cloud of cold atoms, that acts as a 
"classical" dielectric medium of index mo(r). The imaginary part of mo originates in the 
single-atom decay term and is responsible for the diffusive nature of the cloud: it vanishes 
only in the limit of far-detuned incident laser. 



3. Scattered intensity and radiation pressure force 

From the local response /3(r, t) to the external field, many measurable quantities can be 



deriv ed, such as the scattered intensity and the radi ation pressure force (jCourteille et al. 
2010h . The scattered field is (jBachelard et aZ.ll201lh : 

E s (r,t) = -Eo I drV(r0 ^o]r- r1) - 
J fc |r-r'| 

and the far-field scattered intensity, at distance r and direction (9, (f) with respect to the 
z-axis, reads 

J,(r,M) = ceoy^ k(|/3(r)| 2 ) + iV 2 | S (k)| 2 l (3.2) 



where 



(|/3(r)| 2 )^- / dr>(r')|/3(r')| 2 (3.3) 



and s(k) is the structure factor, 

,s(k) = 1 J drp(r)/3(r)e-' ;k r (3.4) 

with k = k (svn.O coscj) : sin sin 0, cos#). The first term of Eq. (|3.2[) corresponds to the 
incoherent and isotropic contribution, and is proportional to N, whereas the second 
term is the superradiant, strongly directional, scattered intensity, and it is proportional 
to N 2 . Hence, the cooperation of an increasing number of atoms to scatter light leads 
to a more coherent and focused emission - superradiance dominates over single-atom 
scattering. 

Another observables of interest is the radiation pressure force, that is the net force 
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that the light exerts on the atoms. It is the sum of two contributions, the first due to 
absorption of photons from the incide nt field and the second one due to their sponta neous 
re-emission along any direction (6, (f>) ( Courteille et al. 201dt Bachelard et aL[2011 ). The 
force on the cloud center-of-mass and along the z-axis reads 



where 



(F z ) = (F?) + (F?), 



E 2 

-27re 01 J-Im( S (ko)), 



-eo 



Ik 2 



dflk cos#|s(k) 



(3.5) 

(3.6) 
(3.7) 



with dtlk = d(j)d9sm9 the elementary solid angle. The absorption force F a is directed 
along the z-axis and it is proportional to the incident intensity, Iq = £qcFq since it 
corresponds to the absorption of laser photons. The emission force F e is directed along 
k, and its z-component is obtained by averaging over the solid angle the product between 
the emission probability of the photons |s(k)| 2 and the projection of their momentum on 
the z-axis, hence the factor cos 9. The integration over the solid angle in Eq. (|3.7|) is 
performed using the identity 



dfi fe cos6>e lk(r - 



-nh(ko\ 



4ni d 
ko dz 



jo(fco|r-r'|), 



(3.8) 



where jo(x) = sin(a;)/a; and j\(x) = sin(a;)/a; 2 — cos(a;)/;r are the zero-th and first-order 
spherical Bessel functions. Using Eq. (|3.8[) . the emission force is written as 



go 
k 2 N 



Re 



drp(r)/3*(r) 



d 



rfrV(rQ eXP , ( ^ o|r - | r,| ^(rQ 



Combining Eqs. 
1 



and 



dz J k \r 
allows to write the total force as 



(3.9) 



N 



drp(r)Re<^ -d^*(r)V 2 



E e 



iko-r 



Eo 



drV(r')- eXP( ^° |r 



r 'l) 5 



fcolr 



/3(r') 



(3.10) 

Since the term in the square bracket is the total electric field E t (r) = i?oe lk °' r + E s (r), 
i.e., the sum of the incident and scattered field given by Eq. (|3.1j) . the force on the center 
of mass appears as the average over the local force on the atoms: 



-dRe(/3*V r £; t )- 



(3.11) 



Hence, w e recover the well-know n expression of the force that a light field exerts on 
an atom ( Gordon fe Ashkin 1980t ). though a crucial difference is that E t here contains 
the self-radiated field of the cloud. Note also that in Eq. (j3?TTJ) , we have extended our 
expression (|3.10|) to every spatial direction. 



4. Mie scattering 

As discussed before, the knowledge of the excitation probability /3(r) is the key to 
understanding the radiation pattern and the forces exerted on the atomic cloud. Despite 
the three-dimensional nature of the problem, an analytical solution exists for simple 
geometries. We here focus on clouds with spherical symmetry p(r), for which the eigen- 
modes of the wave equation have the form u n (r)P n (cos 9), with P n (x) the n-th Legendre 
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polynomial and u n (r) a radial mode that satisfies (jBachelard et al. 2012&1 ): 



<(r) 



k^ml(r) 



n(n + 1) 



u n (r) = 0. 



(4.1) 



Because of the rotational symmetry, no dependence on <j) is present. The excitation field 
is then decomposed as a sum of partial waves: 



0(r,0) = ^(2n + l)^„(r)P Ti (cos^), 

n=0 

whereas the incident wave is decomposed as: 

oo 

e ik °- r - £ i n (2n + l)j n (kor)P„(cos0). 



(4.2) 



(4.3) 



n=0 



Finally, the scattering kernel is expanded in partial waves, using the identity 
exp(zfc |r - r'|) 



oo n 



k \r - 



4 ™E E Y n.m(r)Yl m (f') 



n—0 m— — n 



jn(kQr')h { n\k r) for r > r' ^ 
jn(kor)hn\k r') for r < r 



where f and r' are unit vectors in the directions of r and r', Y n>m the spherical harmonics 
(in particular, Y n fi(9, 4>) = \/ (2n + l)/47rP„(cos#)), and and h„ \z) are the spher- 



ical Bcsscl functions. Inserting Eqs. (14.21) . (|4. 3[) and (|4.4I) into Eq. (|2.4|) and projecting 
on the orthogonal basis of the Legendre polynomial, we obtain the relation 



where 



j n (k Q r) = (28 + i)p n u n (r) + /3„/„(r), (4.5) 
jn(kor')hn\k r) for r > r' 



f n (r)=4m dr'r' 2 p{r')u n {r') 



j n (k r)hn\k r') for r < r' 



(4.6) 



Once specified the form for p(r), the solution for the radial mode u n can be derived, and 
the amplitude of each partial wave j3 n calculated from Eqs. (|4.5I) and (|4.6p. Assuming a 
finite cloud of radius R, such that p(r) = for r > R, it reads (jBachelard et aZ.ll2012M) : 

j n (k R) 



^" (2<5 + i)u n (R) + i\ n hW {k R) 
with 

poo 

X n = 47r / drr 2 p(r)j n (k r)u n (r), 
Jo 

and where we have used that 

/ n (iZ) = 4mh^\k R) / dr'r' 2 p(r')^(r')in(fcor') - i\ n h^\k R). 
Jo 

The cloud structure factor is then easily derived ( Bachelard et al. 201 lh : 

s(k) = — ^2(2n + l)\ n p n P n (cos6), 
as well as the isotropic radiation contribution: 

oo 

(\m\ 2 ) = xJ2( 2n+1 ^\Pn\ 2 > 



(4.7) 
(4.8) 
(4.9) 

(4.10) 
(4.11) 
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Figure 1. Left: Radiation pressure force as a function of the phase-shift $ for homogeneous 
(plain line) and Thomas-Fermi (dashed line) distribution. Simulations realized for a cloud of size 
koR = 10, a detuning 5 = — 50, and varying the number of atoms N. The force is normalized 
by the single atom force Fi — 2iv£o(Eo / ko) 2 / (1 + 4<5 2 ), that is observed in absence of collective 
effects. Right: Radiation pattern 1(6) for homogeneous (thick black line) and Thomas-Fermi 
(blue thin line) distributions. Simulations realized for a cloud for size koR = 10, detuning 
5 = -50 and N = 6100 atoms. 



where 



drr 2 p(r)\u n (r)| 2 



(4.12) 



The problem of scatt ering by a homogeneous spherical dielectrics was initially investi- 
gated by Gustav Mie (Mie 1908). and it was later generalized to homogeneous ellipsoids. 
In the case of a homogeneous spherical atomic cloud of N atoms and radius R, the index 
is m = yl — 3N/(koR) 3 (25 + i), and the radial solutions of the Helmholtz equation is 
u n (r) = ini'mak^r). Using the special properties of the Bessel functions, A„ is explicitly 
calculated: 



A„ = (28 + i)(k a R) 2 [mojn~i(rn k a R)j n (k R) - j n -i(k Q R)j n (m k R)] 
and the complex amplitude of the n-th partial wave is deduced: 



n - 



j„(k R) 



(25 + i)j n (m k R) + i\ n hW(k R) ' 



(4.13) 



(4.14) 



As can be observed in Fig. [I] when the optical density is tuned - here varying the number 
of particles at fixed volume and detuning, the force that pushes the cloud exhibits some 
oscillations. These resonances correspond to the cloud acting as a resonant cavity for the 
light, and are best understood by introducing the phase-shift experienced by the light 
inside the cloud: $ = J(m (0, 0, z) — l)dz. When this phase-shift is a multiple of n, the 
cavity that the cloud forms is at resonance with the light, so a greater amount of light 



is stor ed in the atomic cloud and the radiation pressure force increases (IBachelard et al 
201261 ). 



Yet the spiky structure of the force for homogeneous clouds hides more complex 
resonances. It is known that the sharp index interface between a dielectric and, say, 
vacuum allows for surfac e modes to propagate. These are called whispering gallery 
modes ( Nussenzveid 1992 ). They are characterized by sharp resonances whose number, 
positions and amplitudes strongly depend on the cloud cha racteristics, making them a 
powerful tool to characterize the scatterers ( Oraevskvll2002 ). 
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This exact solution of the scattering problem as an infin ite series of partial waves is 
reminiscent of Mie solution (jMiei Il908l : van de Hulstl 1 1 8 lh , although our approach is 
slightly different: Mie used continuity equations for the tangential field components at 
the dielectric boundary, while we adopted an in tegral formulation of th e problem. The 
two solutions are nevertheless formally the same (jBachelard et aZ.ll2012a> . as far as finite 
clouds are concerned. Let us also remark that our solution is more general, since it holds 
for any radial solution u n (r) of the Hclmholtz equation, i.e., it applies to any spherical 
cloud with arbitrary density. It is thus of particular interest for atomic clouds where the 
traps in general generate inhomogeneous density profiles. 



5. Ultra-cold clouds 

Ultracold clouds of fermionic species have been successfully described by quadratic 
profiles of density, following the pioneering work of Lle wellyn Thomas and Enrico Fermi 



on distributions of electrons (jThomasll 1 927t lFermilH927h . These clouds exhibit a parabolic 



density p{r) = (5JV/2V)[1 - r 2 /R 2 }, with R the radius of the cloud and V = 4irR 3 /3 its 
volume, and consequently a spatially-dependent refraction index: 

mj^r) = m 2 + 7 V, (5.1) 



with m c = y/l — (15/2)iV/(fco-R) 3 (2(5 + i) the index in the core of the sample, and 
7 2 = (l5/2)N/(k 3 l R 5 )(25 + i). Using the substitution u n (r) = r - 3 / 2 w n (x), with x = 
7r 2 /2, one can show that w n (x] satisfies the Coulomb wave equation, well known in 
nuclear physics (jMartinlbOOl) . whose solutions are the Coulomb wave functions w n (x) = 



F n / 2 -i/4 : (—m 2 /'i'y,x). We get the following partial wave expansion for the excitation 
field: 

0(r, 9) = £(2n + ^ P„(cos0). (5.2) 



r 3/2 3-3 I 4 7 ' 2 
n=0 ^ 



Apart from the fact that studying the scattering by atomic clouds with a quadratic 
density rather than an homogeneous one is a much more realistic approximation, the 
former profile also yields different physics. Indeed, as can be observed in Fig.[TJ the spiky 
structure of the whispering gallery modes disappears, and only the regular oscillations 
of the longitudinal cavity modes, such that the longitudinal phase-shift in the cloud is a 
multiple of n, survive. Indeed, because the change in refractive index is much smoother 
in clouds with a Thomas-Fermi distribution than with a homogeneous one, the surface 
modes cannot propagate anymore. 

This analysis is confirmed by the profiles of light intensity in the cloud, see Fig. [21 
One can observe that homogeneous spheres exhibit off-axis surface modes (left picture) 
- this feature was checked for various set of parameters. This surface modes create new 
resonances that compete with the longitudinal ones, causing the spiky structure observed 
in Fig.[T]for the center-of-mass force. These off-axis modes appear much weaker in clouds 
with Thomas- Fermi distribution (right picture in Fig. [2]). 



6. Thermal clouds 

Finally, we discuss the Mie scattering solution in clouds with boundaries extending until 
r = 00. For instance, thermal clouds with a Maxwell-Boltzmann velocity distribution 
have a Gaussian density profile, so the Mie scattering solution that makes use of finite 
boundary conditions cannot be used straightforwardly. Here we propose an alternative 
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Figure 2. Light intensity in the cloud with homogeneous (left) and Thomas-Fermi (right) 
distributions, in the plane y = 0. Simulations realized for a cloud of size koR = 10, a detuning 
5 — —50 and N = 6100 atoms. The intensity is normalized to the laser intensity Iq. 



derivation that holds as well for such distributions with boundaries at infinite. To our 
knowledge, it is the first solution to the Mie scattering problem with infinite boundaries. 

Starting from Eq. (|4.5[) . the amplitude j3 n of the n-th partial wave can be obtained by 
multiplying both terms by 4irr 2 p(r)j n (kor) and integrating over r: 



An 



(25 + i)\ n + B n 



where we have introduced: 



A„ = 4tt / drr p(r)j n (k r), 
Jo 



(6.1) 
(6.2) 

B n = An I drr 2 p(r)j n (k r)f n (r). (6.3) 
Jo 

Using Eq. (|4.6[) . the coefficient B n can be written as 

B n = lQn 2 i J drr 2 p(r)j n (k Q r) x lh^(k r) J dr'r' 2 p(r')u n (r') Jn (k r') 

+ 3n(k Q r) f dr'r ,2 p(r')u n (r')h^(k Q r')} . (6.4) 



Using the identities (|A5[) and (|A 6|) of Appendix \X[ and the identity j n (z)/i„ (z) 
hn (z)j' n (z) — i/z 2 , one can show that 



with 

One deduces that: 



B„ = (26 + i)(-\ n + iA n D n ) 
D n = lim \r 2 \h^(k r)u n (r) - u n (r)h^' (k r)]} 

r—>oo L L J J 



i 



lim 



1 



i(26 + i)D n 26 + i ™ 1 %r2 \ h W {kor)<{r) _ Un ( r )h^ (k Q r) 



(6.5) 
(6.6) 

(6.7) 
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As for the eigenvalues A„ of the scattering problem, they read: 



X n = 47r / drr 2 p(r)j n (k r)u n (r) = (26 + i) lim {r 2 [j n (kor)u' n (r) - u n (r)j' n (k r)]} 



(6.8) 

Thus, the scattering coefficients are determined by the value of the radial modes u n (r) 
and of their first derivative at r — > oo. 



7. Discussion 

We have here discussed the light scattering by a macroscopic atomic cloud, when 
the atoms cooperate to scatter the light superradiantly. The cloud was considered as 
a fluid, i.e., the point-like nature of the microscopic scatterers was neglected, and an 
analytical solution was then derived for spherical geometries, where the excitation field 
inside the cloud is developed as the sum of partial waves. Although we generalized it to 
arbitrary spherical densities, this technique is formally equivalent to Mie scattering, where 
continuity equations are used at the boundaries of the scattering medium. Furthermore, 
our technique allowed for the derivation of a solution to the Mie problem for clouds with 
infinite boundaries. An accurate treatment of the decay of the density profile is crucial 
to understand if some special resonances, such as whispering gallery modes, may exist 
or not. It is also important in the context of atomic clouds or plasmas where, differently 
from solid dielectrics, the densities are usually strongly non- homogeneous. 



Appendix A. 

The function u n (r) is a solution of the differential equation (|4.1j) : 



< + 2^ 
r 



k o m o( r ) - 



n(n + 1) 



u„ = 0. 



(Al) 



where rn^(r) = 1 — An p(r) / k$(26 + i). Defining v n (r) — ru n (r), Eq. (|A 1[) becomes 



k Q m (r) 



n(n + 1) 



v n = 0. 



(A 2) 



For mo(r) = 1, the solutionis q n (r) — rj n (kor). Introducing P(r) — ^mj(r)- n(n+l)/r 2 
and Q(r) = k 2 — n(n + l)/r 2 , we get v'^ + Pv n = and g" + Qq„ = 0, so that 



dr(Q - P)v n q n = q n v' n - v n q' n . 



Since 



(Q - P)v n q n =kl(l- ml)v n q n 



4nrp(r, 2 



k (26 + i) 



r u n (r)j n (k r), 



(A3) 
(A 4) 



we obtain the indefinite integral 



47T 



— / dr'r /2 p(r')u n (r') Jn (k r') = (26 + i)r 2 {j n (k r)u' n (r) - u n (r)j' n (k r)} . (A 5) 



In a similar way we obtain 

47T 



ko 



dr'r' 2 p(r')u n (r')h^(k r') = (26 + i)r 2 [h^ (k r)u' n (r) - u n (r)h^' (k r)) 



(A 6) 
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